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We show how the quark free energy can be determined from a string 
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the relevance of string breaking and avoided crossings. 
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1. Introduction 

What is the maximum temperature where the bulk of QCD thermody¬ 
namics can be described in hadronic terms, with no explicit reference to the 
underlying quarks and gluons ? In this contribution we analyze to what 
extent can one describe the QCD thermodynamics from the hadronic spec¬ 
trum. We thus hope to identify the nature of non-hadronic precursors of the 
crossover to the quark-gluon plasma in terms of low temperature partonic 
expansions m (see [3] for a review). The relation between spectrum and 
partition function provides the thermodynamic free energy F 

Z(T) = X e~ En/T = e~ F,T . (1) 
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The example of a Hydrogen gas in QED illustrates the situation and moti¬ 
vates the discussion regarding the completeness of states which ultimately 
are interacting protons and electrons. At low temperatures, states are 
molecular, i.e. PQ states which behave as compact and structureless. When 
we heat up the system a rovibrational spectrum of states PP|,PP|* emerges 
and eventually the molecule is disociated into atoms H 2 -<->■ 2 H. At higher 
temperatures, atoms are excited H*,H**,... and eventually ionized into a 
plasma of p and e~ in the continuum where the p and e~ constituents be¬ 
comes manifest. Clearly, neutral states (atomic or molecular) are not effec¬ 
tively complete at all temperatures and residual interactions corresponding 
to background non-resonant scattering become relevant [4]. 


2. QCD vs the Hadron Resonance Gas 

In QCD on the lattice F is determined by integrating the trace anomaly, 

e-3P ^ d f P\ 


with respect to the temperature and assuming a reference value, ideally 
P(Tq) = P n (To) with To -C m^. Once this condition is imposed we ex¬ 
pect that at sufficiently low temperatures all the observables are described 
in terms of hadronic degrees of freedom; the quark-gluon underlying con¬ 
stituents are disclosed at very high temperatures. This quark-hadron duality 
at low temperature has been confirmed by the most recent and accurate lat¬ 
tice calculations at finite temperature Eli are confronted with the Hadron 
Resonance Gas (HRG), a rather simple system of non-interacting point-like 
and structureless hadrons 


-4hrg (T) 


jr A I dN{M) I 


d 3 p E(p) - p-V p E(p) 
(2ir) 3 e E ( p ')/ T + T] 


( 3 ) 


Here, the sum is over all hadronic states including spin-isospin and anti¬ 
particle degeneracies, 7 / = =pl for mesons and baryons respectively, E(p) = 
\Jp 1 + M 2 is the energy and the cumulative number of states (degeneracy 
included) N(M ) = ^2 n 0(M — M n ) where M n are the hadron masses. The 
spectrum can be taken either directly from the PDG [7j or the relativized 
quark model (RQM) Refs. [8)0]. Indeed, requiring 


X 2 = 


Alat 

£ 

1=1 L 


-4-Lat {Ti) — -4hrg(T)) 


A*4-Lat (Ti) 


Nda,t ± V 2Adat , 


( 4 ) 


we get T < 175MeV for IVdat = 10 lattice points [5JEJ both for PDG and 
RQM, showing that they can be used to saturate the hadronic spectrum 
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Fig. 1. Left: Total cumulative number N(M) = ^2 n 6(M — M n ) for PDG (full), 
RQM (dotted) and a fit (dashed), Eq. Q with N(M) = A\e M / Ta Q(M t h — M) + 
e Mth / T H0(M — M t h)]. For M t h = 2GeV gives A = 1.5 and Th = 300MeV with 
X' 2 /Y = 0.94 and T max = 185MeV. Right: The single contribution of the deuteron 
3 5i-channel to N(M ) as a function of the invariant NN mass. 


for light quarks provided M < 2.1GeV. This RQM saturation will be as¬ 
sumed also to take place when heavy quarks are involved in which case the 
PDG lacks many of the needed states. This of course suggests to determine 
directly N(M ) from the A data, as illustrated in Fig. |T] (left panel). 

3. Counting hadrons 

The partition function and its comparison with hadronic states through 
the trace anomaly A, displays explicitly the connection between spectrum 
and thermodynamics yielding the statistically significant identification of 
-ZVqcd with A^hrg using PDG or RQM where only qq and qqq states con¬ 
tribute. But, which hadronic states count into the calculation of the ther¬ 
modynamic properties ?. As argued long ago [10J, counting hadronic states 
implicitly averages over some scale, and so states such as the deuteron (made 
of 6 q) generate fluctuations in a smaller nuclear scale. The cumulative num¬ 
ber in a given channel a with bound states M U)a below the threshold Mth 
is in general (see e.g. EQ) 

N a (M) = V 6(M - M n oi ) + - [. 5 a (M ) - 6 a (M th )} , (5) 

z —' 7T 

n 

which becomes N(oo) = ns + [5(oo) — <5(M t h)]/vr = 0 due to Levinson’s 
theorem. This fluctuation is shown in Fig. |T] for the NN channel where 
M t h = 2 Mat, and the deuteron mass is M c ] = 2Mjy — Bd ( B^ = 2.2MeV). 
Thus, the deuteron doesn’t count. This is a general feature of weakly bound 
states, which for the copious new X,Y,Z states, might affect the thermody¬ 
namics if the HRG was blindly identified with the PDG (with all X,Y,Z’s). 
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4. Quark potential, string breaking and avoided crossing 

The linearly growing static energy of two Q and Q sources in the funda¬ 
mental representation of the SU(lV c )-group placed at a distance r is often 
identified with confinement with a string tension a = (0.42GeV) 2 . Including 
the short distance (perturbative) Coulomb-like behaviour yields the vener¬ 
able Cornell potential for the ground state 

A.r\ 

V Q Q(r) = + ( 6 ) 

Actually, the mass of a QQ state, becomes unstable at a critical distance 
r c when a light qq pair is created from the vacuum and the heavy-light 
meson-antimeson BB = ( Qq){Qq ) channel opens 

Mqq ( r c) = Vqq i r c) + mg + rriQ = Mg + Mb , (7) 

where the weak van der Waals like meson exchange interaction in the BB 
sector is neglected 0- We have a diabatic crossing structure which turns 
into an adiabatic avoided crossing when the transition strength Vqq^bb is 

included. In general, one has excited meson states, ^g’^( r ) = A|A + A 

where A ig = A4^ l> — itiq and in Fig. [2] we show the resulting adiabatic 
spectrum based on the RQM |8] when a mixing strength of 50MeV is im¬ 
plemented. 


5. Heavy QQ free energy 


Much of our understanding of the QCD dynamics at finite temperature is 
linked to the Helmholtz free energy defined as the maximum work the system 
can exchange at fixed temperature between two heavy QQ at separation r. 
In the confined phase, the correlator between Polyakov loops becomes [3] 

e -F ave (r,T)/T = (Tr F H(r)Tr F H(0) t ) = ^ e _V 0Q m)(r)/T , (8) 

n,m 


where, before mixing, one has the crossing among the energy levels up to 
qq pair creation. Neglecting the avoided crossing, Eq. (j8]) yields 



1 A rough estimate of r c proceeds as follows. Due to spontaneous chiral symmetry 
breaking the quarks acquire an additive dynamical mass Mo contribution to the total 
mass Mi = Mq + rrii where rm is the current mass. For light qq-mesons M p = 
2Mo + 2 m q whereas for light-heavy (qQ) mesons Mb = 2Mo + rr iq + m,- and thus 
Vgqfrc) = 4Mo + 2m 9 which yields r c ~ 4Mo/cr ~ 1.2— 1.4fm in fair agreement with 
the lattice estimate m- 
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Fig. 2. Left:Potential as a function of distance for the model including heavy-light mesons 
with a charm quark [8], up to A = 1.8 GeV. We include only a g — 50 MeV mixing between 
the fundamental and the excited states. Right: Heavy QQ potential as a function of dis¬ 
tance extracted from the lattice data of Ref. [12] for temperatures T/T c = 0.76, 0.81, 0.90 
and 0.96 using Eq. ©■ The dashed (red) line is the result for the QQ Cornell potential 
of Eq. (O with a = (0.4 GeV) 2 . 


where A n = AjA = A^q. We show in Fig. [3] the heavy QQ free energy 
obtained with Eq. d9]) by using the spectrum of heavy-light mesons with a 
charm (bottom) quark, and no strangeness, obtained with the Isgur model 
of Ref. |8l up to A = 3.19 GeV. We have considered a = (0.4 GeV) 2 and 
a = 7t/ 16 (no fit). In Fig. [3]we also show the Polyakov loop computed as 

L(T ) = lirn e -Fwe(r,T)/(2T) 1 V' e -A„/T _ ( 10 ) 

r—>oo 2 J 

n 

The agreement with lattice data is quite good for T < 0.8 T c . Thus the 
spectrum both for the PDG and the RQM saturates the sum rules at these 
temperatures @ whereas the finite temperature does not resolve avoided 
crossings making them effectively diabatic passages. 


6. Conclusions 

While much of the effort has been directed towards a description of the 
hadron-to-quark-gluon crossover, at present there is no understanding on the 
mechanism, perhaps because confinement itself is not really well understood. 
Our analysis of the free energy suggests that lowest temperatures where 

2 By multiplying L(T) by a factor e c ^ T with C = —40MeV, we get good agreement 
with the lowest temperature lattice point. This kind of ambiguity comes from renor¬ 
malization effects. Like in | 12 | we assume Ekve(r, T) ~ LqqR) for rT 1. The 
recent evaluation of Ref. ns. unlike m is based on a different renormalization con¬ 
dition and Nf = 2 + 1 for the free energy and the role of avoiding crossing will be 
analyzed elsewhere. 
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Fig. 3. Left: Color averaged heavy QQ free energy as a function of distance. We 
show as dots the lattice data for Nf = 2, T = 0.76 T c taken from Ref. j!2j . and 
as continuous lines the result by using Eq. Q with the spectrum of heavy-light 
mesons with a charm quark (red) and bottom quark (blue), with no strangeness, 
obtained with the Isgur model of Ref. [8], as well the spectrum obtained from the 
PDG (green). Right: Renormalized Polyakov loop as a function of temperature. 
We show as dots the lattice data for T = 0.76T C taken from Ref. 02, and as 
continuous lines the result by using Eq. m- 


purely hadronic states account for QCD observables is rather large T ~ 
170 — 180MeV and close to the critical temperature. 
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